Greenberger-Horne-Zeilinger States in Quantum Dot Molecule 
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We present a microscopic theory of a lateral quantum dot molecule in a radial magnetic field 
with a Greenberger- Horne- Zeilinger (GHZ) maximally entangled three particle groundstate. The 
quantum dot molecule consists of three quantum dots with one electron spin each forming a central 
equilateral triangle. The anti-ferromagnetic spin-spin interaction is changed to the ferromagnetic 
interaction by additional doubly occupied quantum dots, one dot near each side of a triangle. The 
magnetic field is provided by micro-magnets. The interaction among the electrons is described within 
an extended Hubbard Hamiltonian and electronic states are obtained using configuration interaction 
approach. The set of parameters is established for which the ground state of the molecule in a radial 
magnetic field is well approximated by a GHZ state. 

PACS numbers: 03.65.Aa, 03.65.Ud, 03.67.Bg, 71.10.Fd, 73.21.La 



I. INTRODUCTION 

In a recent work Roethlisberger et al., Ref. HlI, proposed 
a scheme to generate maximally entangled Greenberger- 
Horne- Zeilinger— (GHZ) state as a ground state of a 
three spin system. The spins were arranged in a triangu- 
lar geometry and assumed to interact ferromagnetically. 
With ferromagnetic interaction the degenerate ground 
state is maximally spin polarized, with two | ttt> and 
I 4,4,4, > configurations. When radial in- plane magnetic 
field is applied the two spin polarized configurations form 
the two GHZ states \GHZ^ >= -^{\ ttt> ±1 iii>)- 
Such GHZ ground states of spin molecules could be 
used as long lived sources of entanglement^ii. One can 
envisage using lateral triple quantum dot molecules^"— 
with one electron on each dot to realize a three spin 
system. However, the spin-spin interaction in quantum 
dot molecules is necessarily anti-ferromagnetic and such 
a simple quantum dot molecule is not possible. Here 
we propose more complex molecule which uses auxil- 
iary quantum dots to effectively change the spin-spin 
interaction to ferromagnetic and generate the maximally 
entangled GHZ state as its ground state. 



of each pair of localized spins is modified by connecting 
each pair of spins to an auxiliary quantum dot with two 
electrons. In Ref. [l3 it was shown that the effective 
interaction between two electrons in two quantum dots 
i.e., one electron in each of the two quantum dots, can 
be changed from anti- ferromagnetic to ferromagnetic by 
connecting them to a third (auxiliary) doubly occupied 
quantum dot and applying a bias to the auxiliary quan- 
tum dot. Thus it was shown that there exists a set of 
parameters for which such a triple quantum dot molecule 
with four electrons can be thought of as effectively two 
localized electron spins interacting ferromagnetically. 
Following this idea we propose here a quantum dot 
molecule consisting of three dots with one electron each 
and three auxiliary dots with two electrons each which ef- 
fectively realizes the ferromagnetically coupled three spin 
cluster. Using Hubbard model and exact diagonalization 
techniques we show that the ground state of such a 
molecule in a radial magnetic field is indeed a GHZ state. 
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Theoretical proposals to realize GHZ states in solid 
state systems include, e.g., spin systems^, excitons in 
coupled dotgSi, two-level atoms in a non-resonant cavity^ 
and superconducting flux qubits^*^. The GHZ states have 
been experimentally realized using photonsii, in atomic 
systems using three Rydberg atomsi^ and very recently 
maximally entangled GHZ states have been realized in 
solid state system using superconducting qubits^'^'^'*. 
We refer the reader to Refs. [islfl^ for full initialization, 
characterization and read out of GHZ states. 

In this paper, we propose a theoretical model of six 
coupled quantum dots-i^ '^^'^^ with interacting electrons 
to generate the GHZ state as a ground state. The model 
consists of three equally spaced quantum dots with one 
electron spin each. The anti-ferromagnetic interaction 



The structure of the paper is organized as follows. In 
Section [IT] we examine the proposed setup for generating 
highly entangled states in a three spin system interact- 
ing ferromagnetically with each other and coupled to 
a radial in- plane magnetic field. We clarify the role 
of applied magnetic field in generating the GHZ states 
using the degenerate perturbation theory. In Section IIIII 
we propose a quantum dot molecule consisting of six 
quantum dots with nine electrons and micro-magnets 
to generate radial in- plane magnetic field creating the 
maximally entangled GHZ state as its ground state. We 
describe the system, its Hamiltonian and ground state 
phase diagram. The last section IIVI contains summary 
and conclusions. 
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II. GENERATING GHZ STATE IN A THREE 
SPIN SYSTEM 

We begin with an analysis of the three spin system 
to create the maximahy entangled states as outlined in 
Ref. [l|. The three spins | placed at the corners of an 
equilateral triangle lying in the xy plane as shown in left 
panel of Fig. [1] are described by an isotropic Heisenberg 
Hamiltonian, 



eigenstate, S the total spin and Sj the z-component of 
total spin of the system. 

The ground state 5 = 3/2 quadruplet consists of the 
following four eigenstates with eigenvalue — |J, 

1,10 = I ^tt> (2) 



(S4 — Si^ 



(1) 



(3) 



where it is assumed that the exchange coupling J is 
ferromagnetic (J > 0) and Si — ^(crf , erf , erf ) is a vec- 



tor consisting of Pauli matrices, af 



1 

1 



^ and ''■f — ^ Q ^ij ^*^1^^^S "-"^ ^^^"^^ ^ (=1,2,3). 

Here h is reduced Planck constant. 
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FIG. 1: (Color online) Schematic showing two views of three 
localized electron spins Si (i=l,2,3) arranged in a triangular 
geometry and interacting ferromagnetically with each other. 
The two bar magnets at the triangle centroid are placed on 
top of each other generating a radial magnetic field Bi. Figure 
on the left shows the top view of such a system and the right 
panel exhibits the lateral view. The magnetic field lines are 
also shown. 

We write the Hamiltonian, given in Eq. [TJ in the basis 
of tensor product of states at each site. For example, the 
Pauli matrices are written in such a local basis. Thus 
the Hamiltonian in the basis set {| ttt>, I tti>7 1 tlt> 
, I til>, I itt>, I iti>, I lit> and| m> } is a 8 X 8 
matrix. 

Exact diagonalization of the Hamiltonian matrix 
in the basis of 8 states yields two quadruplet levels 
corresponding to total spin 5 = 3/2 and total spin 
S = 1/2. We write the eigenstates of the Hamiltonian 
in Ec^. [1] as |a, 8,8^ > where a is the index of the 



3,|0 = ^[un>+int>+im>] (4) 



'2' 2 



^ = -^[im>+un>+u;t>] (5) 



The four states correspond to the two distinct classesi^ 
of highly entangled states of a tripartite qubit system. 
As will be shown below, the mixed |1> and |2> states 
form the two states of the Greenberger-Horne-Zeilinger 
class [|GHZ± >= ^(1 ttt> ±1 iU> )] while |3> and 
|4> belong to the W clas&i^. 



The excited quadruplet with total spin 5=1/2 and 
genvalue 
states^^. 



eigenvalue comprises of the following four chiral 



1 1 

^'2' 2 
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[I ;n> +e"^ I m> +e"^ I m> ] (6) 



6,i,i\ = 4l 



2' 2 



1 1 
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;tt> +e ^ I tlt> +c 



m> 



(7) 



V3 



71' 



[I n;> +eT^i in> +eT^i ut> 



(8) 



(9) 

These chiral states are the eigenstates of chirality 
operator x = Si • (S2 x S3) as described in Ref. H^. The 
states have the spin ( f and 4, ) currents going in two 
different ( and O ) directions. 

For simplicity we shall shift the overall energy scale 
by I J with the energy of the ground state quadruplet 
as the reference energy. In order to obtain one of the 
two GHZ states as a ground state of the system we need 
to apply a perturbation which will mix the states |1> 
and |2> while separating them from the remaining two 
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W states, |3> and |4>, thus splitting the degenerate 
eigenspace of the lower quadruplet. As was shown in 
Ref.lH this can be accomplished by applying an in- plane 
radial magnetic field. 

We now apply the radial magnetic field B as shown in 
Fig- HI The Hamiltonian Hb = go/^B X^i -^i ' describing 
the coupling of the spins to an externally applied in- plane 
radial magnetic field generated using the two bar magnets 
placed with similar poles2i facing each other as shown in 
the right panel of Fig. [TJ reads 



Hb =b 



(T9 — 

2 



2 ^ 



V3 



(10) 

where b — — - is an effective magnetic field with g^, 
an effective electron g- factor in the plane of a quantum 
dot, /zb the Bohr magneton and B the strength of the 
radial magnetic field at the quantum dot position. The 
effective magnetic field 'b' is simply related to Zeeman 
energy splitting 2b — AE^. for a given external magnetic 
field 'B'. And Sj is a vector denoting the spin at quantum 
dot "i" with Pauli matrices as components along x and 
y directions. 



H 




(13) 



The effective Hamiltonian is that of a two level system 
with tunneling T — proportional to the third power 
of the external magnetic field 'b' . After diagonalization 
of the matrix in Eq. [T3] we obtain as the ground state 

the GHZ state \GHZ+ >= ^(| ttt> +1 IU> with 

energy — fjr, separated from the second GHZ state by 
2T = ^ as also obtained in Ref. U 

The energy splitting of the two GHZ states can be 
also written in terms of exchange coupling J and the 
Zeeman energy as 2T = j{AE.^){^)^. For GaAs 
AEz=20 fieV in the field of B=l Tes'la. The exchange 
coupling measured for coupled lateral GaAs quantum 
dots is comparable to Zeeman splitting as J^20 /ieV— . 
Hence in lateral devices the two GHZ states can be 
separated on the order of Zeeman energy, i.e., several ^ 
eV. The full discussion of physical parameters in gated 
and self-assembled quantum dot is deferred to the end 
of Section Uni 



The complete Hamiltonian of the three spin system in 
an external magnetic field is given as 



H = H 



H 



Hp 



(11) 



The physical meaning of the radial magnetic field is 
best described by its action on the spin polarized state: 



In the next section we discuss how spins of electrons 
in a quantum dot molecule can be used to realize the 
GHZ states in solid state quantum dot molecules. 



III. QUANTUM DOT MOLECULE 



hbI ttt>= ^[i m> +e-'^\ m> +e-'^\ m>] 

(12) 

Thus the radial in- plane magnetic field simultaneously 
flips the spins and adds a phase factor in such a way as 
to generate a chiral state |6 > with S = 1/2,8^ = +1/2 
and momentum fc = 27r/3. In a similar way Hb couples 
state |6 > (St +1/2, k = 27r/3) with state |7 > ( S^ = 
— l/2,fc = — 27r/3 ) and state |7 > with spin polarized 
state |2 > with (S^ = -3/2, fc = 0). The transition 
from state |1 > to state |2 > is a third order process in b. 
The final effective Hamiltonian for the two spin polarized 
states (|1> and |2>) and the two chirality states (|6> and 
|7>) reads, 

/ b%/3 \ 
OOP hV3 
hVS |J 2b 
\ bV3 2b |J / 

With the radial magnetic field (second term in Eq. [TT|) 
treated as a perturbation, the degenerate perturbation 
theorySS, leads to an effective 2x2 Hamiltonian matrix 
^/fe//-j ^Yie subspace of spin polarized S=| states. 



A lateral triple quantum dot molecule 
(TQDM)^"Iiiii2i-^ with one electron per dot is the 
simplest realization of the three spin system discussed 
in the previous section. Assuming a single orbital per 
dot and arbitrary occupation, as shown in Refs. 00 the 
quantum dot molecule can be described by an extended 
Hubbard Hamiltonian, 

N N N 

^QDM ^ J2 e^lMr - J2 (4rCjr +h.C.)+Vj2 n-Hi-r 

i=l i,j = l i=l 

T r r 

N 

+ J2 "^^iP^Pi (14) 

where N is the number quantum dots ( N= 3 for 
TQDM), Cjtj. (Cit) is the creation (annihilation) operator 
of an electron with spin r(=t,-l.) at i-th QD , niT-(= 
cj^CiT-) is the spin- dependent electron occupation 
number, Pii—J^T^iT-) i^ total occupation number 
at site "i", ei is the on- site energy, ty is the inter- dot 
tunneling matrix element between nearest neighboring 
dots "i" and "j", Ui is the intra- dot and Vy the inter- 
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dot Coulomb repulsion between the nearest neighboring 
dots and h.c. means hermitian conjugate. 

At half-filling the Hubbard Hamiltonian reduces to 
an isotropic Heisenberg Hamiltonian but with anti- 
ferromagnetic exchange coupling^iZ,, ~ 4t^/U, and the 
ground state (GS) is in the subspace of total spin S=^. 
Thus it is not possible to tune the GS of the system 
to a maximally spin polarized subspace of where 
the GHZ states reside. However, it was shown^iii that 
triple quantum dot molecule with four electrons, which 
corresponds to two holes, has a spin polarized ground 
state. When one of the dots is biased, it contains two 
spin singlet electrons while the two remaining elec- 
trons, localized on the two remaining dots, are coupled 
ferromagnetically. Biasing the dot further effectively 
decouples this dot and leaves the remaining two dots 
with one electron each in a spin singlet state, i.e., 
coupled anti-ferromagnetically. 

Hence it is in principle possible to design a half- filled 
TQDM combined with auxiliary doubly occupied quan- 
tum dots in order to generate an effective ferromagnetic 
coupling of spins in the TQDM. If this is accomplished 
then we can realize the GHZ states in such a quantum 
dot molecule by applying a radial in- plane magnetic 
field as explained in the previous section. 

In what follows we propose two design structures, 
namely a six and a four quantum dot molecule as shown 
on the left and right side of Fig. [2] respectively. While 
studied both, we discuss in detail the six quantum 
dot molecule because it might be difficult to fabricate 
a triangular quantum dot structure with a central 
auxiliary dot in the middle. 




FIG. 2: (Color online) Left: Schematic of a six dot molecule 
with dots labelled (1,2,3) forming the central triple quantum 
dot molecule and dots (4,5,6) are the auxiliary dots. Right: 
A four dot molecule structure with dot 4 being an additional 
dot introduced to the triple quantum dot molecule (1,2,3). 



We form a six quantum dot molecule (SQDM) by 
bringing three doubly occupied dots labelled (4,5,6) 
close to the central^l triangle (1,2,3) of singly occupied 
dots, for a total of nine electrons. The Hamiltonian for 
SQDM, "H^QD^ gi^gjj Eq.IHfor N=6. We consider 
tij=t and Vij=V for central triangle {i,j} = {1,2,3} and 
tij=t' and Vij=V' for all other dots as shown in left 
panel of Fig. [5] 

We create all 220 No = 9 electron configurations 
J^jg.c^j^|0 >, construct the Hamiltonian matrix and 
diagonalize it numerically to obtain the energy spectrum 
and eigenvectors. In the absence of the magnetic field 
the z- component of the total spin is a good quantum 
number giving the Hamiltonian matrix in a block 
diagonal form. The Hilbert space dimensions of the 
blocks in the subspace of z- component of total spin 
St=±i and St=±| are 90 and 20 respectively. Each 
block is numerically diagonalized so as to obtain the GS 
of the system as a function of the Hubbard parameters. 

In our calculations we assume ci = €2 = — 
i.e., the dots (1,2,3) are on resonance. The energy 
£4 = £5 = £6 = e of auxiliary dots is varied by applied 
gate voltage. We consider the regime U ^ t,V with 
constant value of U=2.0, V=0.1 as considered earlier—. 
We vary inter- dot Coulomb repulsion and tunneling 
between central and auxiliary dots as < V < V and 
< t' < t. All the parameters are in the unit of effective 
Rydberg defined by Ry — m*e'*/2e^/i^ where m* is the 
electron effective mass, e the electron charge, e the 
dielectric constant of a material. For example in case of 
GaAs the effective Rydberg is estimated to be about 6 
meV. 
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FIG. 3: (Color online) Dominant charge configuration in the 
subspace of total spin S=| and z- component of total spin 
St = | for a six dot molecule with nine electrons. 
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The GS of the system is a hnear superposition of the 
basis used to construct the Hamihonian matrix in that 
given subspace. For certain parameters the GS contains 
the spin polarized configuration with spins locahzed on 
the three dots. This is illustrated in Fig. [3] which shows 
the dominant configuration, with probability 0.9289 
for the GS in the subspace of total spin S=| and z- 
component of total spin S^—^ obtained for parameters 
e = -0.05, U=2.0, V=0.1, t=0.05, V'=0.094, t'=0.05. 
It is clearly seen that the three electron spins are 
localized on the triangle of dots (1,2,3) forming the 
ferromagnetically coupled three spin system and the 
auxiliary dots (4, 5, 6) are doubly occupied, driving the 
GS of the quantum dot molecule to be in the maximally 
spin polarized subspace. 



of GS in maximally spin polarized subspace increases. 
It is because increasing the tunneling, t, delocalizes the 
electrons in the central TQD molecule which decreases 
the contribution of the spin polarized configuration, 
Eq. [15] But in this scenario even a small value of 
tunneling t' between the central and auxiliary dots can 
drive the GS in the maximally spin polarized subspace 
and so the area comprising of GS in total spin S=| 
increases. 

It is also observed that the maximum value of overlap 
probability lies along the border line dividing the regions 
where the GS is in the subspace of total spin S=^ 
(zero probability) and that of total spin S=| (finite 
probability) . 



We define a state corresponding to the dominant 
charge configuration in the subspace of z- component of 
total spin 8^=1 as. 



|US I(tlt2t3)(t4;4t5;5t6;6) >, 



(15) 



Choosing an appropriate set of parameters can maxi- 
mize the overlap probability and prepare the GS of the 
system with spin polarized charge configuration given in 
Eq. We thus define a GHZ state for the SQDM sys- 
tem as. 



with similar definition for the dominant charge con- 
figuration in the subspace of z- component of total spin 



IDS I(;i;2;3)(t4;4t5;5t6;6) >, 



(16) 



We now determine the set of parameters for which 
the SQDM describes the three spin system well, i.e. we 
calculate the overlap probability, at zero magnetic field 
(B = 0), 



Pb=o = I < GS|US > 



(17) 



of the numerically obtained GS with state |US) as a 
function of parameters. The non-zero overlap probability 
as a function of parameters t,t',V' and e for fixed U,V, 
determines the phase diagram, region in parameter space 
where the three electrons are localized on the three dots 
and interact ferromagnetically. The result of numerical 
calculation for U=2.0 and V=0.1 is shown in Fig.U] 

We observe that for a fixed value of t, if we vary e 
the region of GS in maximally spin polarized subspace 
gets reduced but the maximum value of the overlap 
probability increases marginally. This happens because 
as we lower the on- site energies of the auxiliary dots, 
the three central dots become isolated from the auxiliary 
system and behave like half- filled TQDM. And so a 
higher value of t' is required to have the ground state 
in maximally spin polarized subspace. Hence the region 
belonging to the total spin S= 



:| decreases. 



It is also found that for a fixed e, varying t decreases 
the maximum value of overlap probability but the region 



\gnz^ >= ^[|us > ±|DS >]. 

v2 



(18) 



In order to form one of the states given in Eq. [TH] as 
a GS of the SQDM we apply a radial in- plane magnetic 
field as shown in Fig[SJ The Hamiltonian describing this 
interaction of the electron spins with external magnetic 
field is, — gc/^B X]f=i -^i ■ ^n-d is given as 



H^=b 



V3 



2 



2 ^ 2 
b 



+ b 
73 



2"3 



V3 y 

2 ^ 



+ b 






(1 



The spin operators (cr'^, a^) are written in second quan- 
tized form as 



and 



(20) 



(21) 



We write the complete Hamiltonian for SQDM and an 
applied external magnetic field as. 



(22) 



and the Hamiltonian matrix is constructed using 
the same basis as used to create Ti,^'^^ . But since 
the magnetic field doesn't preserve the spin rotational 
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FIG. 4: (Color online) Phase diagram of ground state (GS) overlap probability with a state of dominant charge configuration 
in the subspace of S=| for a six quantum dot molecule is shown as a function of t'/t and V'/V. The individual view graphs 
are for different values of hopping element, t and the on- site energy, e, at dots (4,5,6). The lower area in each graph (green in 
color) denotes the GS is in subspace of total spin S=^ and its overlap is zero. The number in each graph corresponds to the 
maximum value of the overlap probability at the point given in the bracket. 



symmetry the Hamiltonian matrix in Eq. [5^ is not in a 
block diagonal form. The dimension of the full Hilbert 
space is 220. We diagonalize the Hamiltonian matrix for 
three different values of e and tuneling matrix element, t. 
These are chosen such that in the absence of an external 
magnetic field the area covered by the subspace of total 
spin S=| in each view graph of Fig. |4] is i) minimum 
[e=-0.10,t=0.04)], ii) intermediate [e=-0.05,t=0.05] and 
iii) maximum [e=0.00,t=0.06]. And for each value of e 
and t, we diagonalize the Hamiltonian matrix for the 
parameters ( ^ , 7 ) where the overlap probability, Eq. [T71 
is maximum. 

At these points the low energy spectrum of extended 
Hubbard Hamiltonian resembles the energy spectrum 
of an isotropic Heisenberg Hamiltonian with four- fold 
degenerate ground and excited state as discussed in the 
previous section. For zero magnetic field the energy 
separation between these two degenerate state is ^ . We 
use this to find the values of ferromagnetic exchange, 
J, in our Hubbard model for parameters (t,t',V,V',e) 
mentioned above with fixed U,V and we use it as our 
energy scale. 



For each set of parameters, exact diagonalization of 
Eq. [21] yields the two GHZ states as the ground and first 
excited states of our system. In Section [H] for Heisenberg 
system, we have seen that the energy separation between 
the two GHZ states for finite magnetic field is ~ ^^^||- 
where b is the effective strength of magnetic field. Apart 
from Heisenberg model, we calculate the energy splitting 
between the two GHZ states for our Hubbard model as 
a function of ^ for three different values of e, t and 
corresponding J. The result is shown in Fig. HI We 
observe that for smaller values of j, Heisenberg and 
Hubbard model behaves similarly but for larger values 
the agreement is worse. But it is also seen (inset of 
Fig. [S]) that as we keep on decreasing e and t the energy 
gap obtained from the Hubbard model starts approach- 
ing the Heisenberg model. This is expected since the 
value of maximum overlap probability is smallest for 
(e=0.00,t=0.06) and largest for (e=-0.10,t=0.04). 

Using the same set of Hubbard parameters, we also 
evaluate the maximum overlap probability of the ground 
state with one of the GHZ state for the SQDM as given 
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FIG. 5; (Color online) Schematic depicting top bar magnets 
(rectangular box shown in blue) with like poles facing each 
other are placed at the center of the system which generates 
a radial in- plane magnetic field. 




FIG. 6: (Color online) Energy splitting between the two GHZ 
states is shown as a function of ^ for three different values 
of e, t and J. The stars represent the results as obtained for 
three spin system, Heisenberg model. 

in Eq. [18] for finite magnetic field and denote it as 

Pb^o = I < GS\gHZ+ > |2 (23) 

The results are shown in Fig. [7] as a function of e, 
t and J. We find that as the strength of the effective 
magnetic field, b, approaches the value of J the GS of 
the Hubbard model deviates more and more from the 



1 



FIG. 7: (Color online) Maximum overlap probability of the 
ground state with QT-LZ^ is shown as a function of 2^ for 
three different values of e, t and J. 



GHZ state for SQDM which results in the decrease of 
the overlap probability. 



We now turn to the discussion of physical parame- 
ters and structures needed for the realization of GHZ 
generator. Typical parameters leading to an effective 
ferromagnetic coupling of three spins localized on a 
triangle, shown in Fig. 21 involve U=2, V=0.1, t=0.05, 
t'<t, V'<V and a bias e of the order of tunneling, t. 
The values currently available for lateral gated quantum 
dots on GaAs are in the required parameter ranges of 
U--2meV, V~0.1meV, and t-0.05meV. By building 
quantum dot networks using individually gated self- 
assembled quantum dots on nanotemplates— one can 
envisage reaching values of parameters for self-assembled 
quantum dots U~20meV, V^lOmeV, t^^lOmeV which 
should lead exchange coupling reaching J^40meV, 
exceeding room temperature. 



Finally, we discuss the lifetime of a maximally entan- 
gled GHZ state. Since in the proposed scheme GHZ 
state is the ground state, it is expected to only suffer 
decoherence and no population decay. If the quantum 
dot molecule is realized in GaAs, nuclear spins are 
expected to be the major source of decoherence. Since 
the problem of decoherence of the GHZ state is common 
with a single electron spin, recently developed coherent 
control of nuclear spins^ might be expected to be 
applicable in extending coherence of the GHZ state. 
This problem will be investigated in the future. 
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IV. SUMMARY & CONCLUSION 

In summary, we designed and analyzed theoretically 
a lateral quantum dot molecule combined with a 
micro-magnet generating a maximally entangled three 
particle GHZ groundstate. The quantum dot molecule 
consists of three quantum dots with one electron spin 
each forming a central equilateral triangle. The anti- 
ferromagnetic spin-spin interaction is changed to the 
ferromagnetic interaction by additional doubly occupied 
quantum dots, one dot near each side of a triangle. 
Exact diagonalization studies of the Hubbard model of 
the molecule determine the phase diagram in parameter 



space and a set of parameters is established for which 
the ground state of the molecule in a radial magnetic 
field is well approximated by a GHZ state. 
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